Strongly renormalized quasi-two-dimensional electron gas in a heterostructure with 

correlation effects 
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We present aspects of strong electron correlation in a model of a BI/M/BI heterostructure where 
BI is a band insulator and M can be tuned from a metal to a Mott insulator by varying the on-site 
repulsion. An effective one-dimensional Schrodinger equation for the low-energy part is derived in 
the framework of the Kotliar-Ruckenstein slave-boson mean-field theory and solved self-consistently 
for a wide range of the on-site interaction strength and of the width of the sandwiched material M. In 
the strongly interacting limit coherent quasiparticles responsible for metallic behavior are confined 
to a relatively narrow region near the interfaces where they form a quasi-two-dimensional electron 
gas. Due to the lack of spatial homogeneity these quasiparticles intrinsically obtain a nontrivial 
dependence on the transverse momentum which for example enters the low-energy behavior of the 
optical conductivity. 

PACS numbers: 71.27. +a Strongly correlated electron systems; heavy fermions, 73.20.-r Electron states at 
surfaces and interfaces, 73.40.-c Electronic transport in interface structures, 78. 20. -e Optical properties of 
bulk materials and thin films 



I. INTRODUCTION 

Strongly correlated electron systems are strong can- 
didates for the next generation functional materials for 
various purposes. This bright perspective is motivated by 
the richness of their phase diagrams involving a variety 
of ordered phases in charge, spin and orbital degrees of 
freedom, metal to insulator transitions, unconventional 
and high-Tc superconductivity as well as various exotic 
liquid phases i The electronic and magnetic phases show 
a high sensitivity to external conditions such as tempera- 
ture, pressure, doping and application of electric or mag- 
netic fields. In artificially built hcterostructures involving 
different strongly correlated electron systems the condi- 
tions may change locally leading to new physics at the 
interface as compared to the bulk material. These sys- 
tems are likely to have remarkable properties which may 
be tuned by the application of electric fields 2 - or by mod- 
ulation doping near the interfaces^ Taking advantage of 
new mechanisms at interfaces of correlated electron sys- 
tems may allow to design junctions for potential device 
applications, e.g., a colossal electroresistance device^ 

One class of recently investigated heterostructures is 
made of complex oxides which arc insulating in the bulk 
systems In some cases, when brought together, a 
metallic quasi-two-dimensional electron gas at the inter- 
face is formed. Pioneering experimental work by Ohtomo 
and coworkers are performed at LaTi0 3 /SrTi0 3 (Ref.0) 
and at LaAlC>3 / SrTiC>3 (Ref. 0) interfaces and superlat- 
tices. In both cases the bulk materials are insulating but 
in the epitaxial hcterostructures the transport properties 
can change dramatically. Thiel et al£ reported on a very 
large electric field response of the quasi-two-dimensional 
electron gas at the LaAlOs/SrTiOa interface which is of 
potential interest for device applications. 

Our work is mainly inspired by the experiments on the 



LaTiOa/SrTiOa interfaced^ As pointed out by Okamoto 
and Millis^ the atomically precise fabrication, the near 
lattice match (a w 3.9 A) and the chemical similarity of 
the two components offer a good starting point for the- 
oretical studies and make it possible to study the influ- 
ence of electronic reconstruction alone: how does the elec- 
tronic phase at the interface differ from that in the bulk 
system? The pure La-compound is a Mott insulator with 
a Ti-Sd 1 configuration whereas the pure Sr-compound is 
a band insulator (Ti-3e?° configuration). When brought 
into contact the mutual doping of band and Mott insula- 
tor leads to a quasi-two-dimensional electron gas confined 
in an approximately 2 nm thick region at the interfaced 
Mainly for single-La-layer superlattices, where the effect 
of strong electron correlation is believed to be minor, 
density-functional theory calculations were successfully 
performe d 10 ! 11 ' 12 : 13 to primarily study lattice-relaxation 
effects and correlation-driven spin and charge ordering. 
In order to include additional effects of electron corre- 
lation Okamoto and Millis proposed a generalized Hub- 
bard model for this system and combined insights from 
the Hartree approximation^*!! and dynamical mean-field 
theory (DMFT)i 15 ! 16 Focus was put on the electronic 
charge profile which was determined in the experiment 
by electron energy loss measurements and on possible 
ferromagnetic and antifcrromagnctic ordering. 

The above discussed systems built an oxide analogy to 
more conventional semiconductor heterostructures. But 
strong on-site interaction in combination with a spatially 
nonuniform setup is expected to lead to differences as 
compared to both the strongly correlated bulk systems 
or heterostructures where electron correlations play a mi- 
nor role. In order to give a systematic discussion on the 
influence of the on-site interaction we study a simpli- 
fied model for hcterostructures which are characterized 
by a BI/M/BI stacking. Here BI is a band-insulator 
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such as SrTiC>3. By varying the on-site repulsion wc 
can continuously change the properties of the sandwiched 
material M tuning it from a metal to a Mott insulator 
such as LaTiO,3. The model is based on a generalized 
single-band Hubbard model introduced in Ref. [TH. So 
far, it was studied within the Hartree^ Thomas-Fermi 3 -, 
and DMFT approximation! 3 ^ While the Hartree and 
Thomas-Fermi approximation are expected to give rea- 
sonable results for the charge distribution, i.e. the self- 
consistent screening, these approaches can not take into 
account the renormalization of the quasiparticle weight 
and some aspects of correlation are thus treated poorly. 
DMFT on the other hand is a powerful tool to describe 
these aspects but especially for spatially nonuniform sys- 
tems a full treatment becomes numerically very involved. 
Here we report on an alternative approach to investigate 
the quasiparticles in such systems, namely, by means of 
the slave-boson mean-field (SBMF) approximation. The 
Kotliar-Ruckenstein (KR) slave-boson theory was origi- 
nally formulated for the single-band Hubbard models It 
reproduces the Gutzwiller approximation^ on the mean- 
field level and is thus able to catch important aspects of 
strong correlation such as the suppression of the quasi- 
particle weight and at half filling the Brinkmann-Ricc 
transition^ from a paramagnetic metal to a param- 
agnetic insulator. Generalizations to a manifest spin- 
rotation invariant slave-boson representation were made 
by Li et al^ and to a unified spin- and charge-rotation 
invariant formulation by Fresard and Wdlflo 2 -^ leading to 
a more profound understanding of the KR slave-boson 
theory. Applications to the orbital degenerated Hubbard 
model were performed to study (anti)ferromagnetis m 22 i 23 
and the orbital-selective Mott transition^ 

The paper is organized as follows. In Sec. [IT] we intro- 
duce the model and apply the SBMF approximation. In 
Sec. Qll] we discuss the T = properties of the param- 
agnetic phase with focus on the coherent quasiparticles. 
Due to the lack of spatial homogeneity the quasiparti- 
cles obtain a non trivial momentum dependence which is 
studied by calculating the renormalization of the quasi- 
particle dispersion and the optical conductivity. In the 
last section a summary and conclusions are given. 



II. MODEL AND FORMALISM 

A. Model of a BI/M/BI heterostructure 

We study the model of a [001] BI/M/BI heterostruc- 
ture where BI is a band insulator and M a material which 
can be tuned from a metal to a Mott insulator by varying 
the on-site repulsion. As compared to the experimentally 
realized LaTiOs/SrTiOa heterostructures the basic sim- 
plifications of the model are the disregard of the orbital 
degeneracy and the assumption of an unrelaxed perfect 
lattice match. A sketch of the studied system is shown in 
Fig. HJ The heterostructure is defined by the electrostatic 
potential of N ion-layers with simple positively charged 
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FIG. 1: Sketch of the heterostructure studied in this paper. 
BI is a band insulator and M is a material which can be tuned 
from a paramagnetic metal to a paramagnetic Mott insulator 
by varying the on-site repulsion. The heterostructure is de- 
fined by N positively charged counter ions sitting in between 
the electronic sites. 



ions sitting in between the electronic sites of a simple 
cubic lattice 
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These counter ions simulate the electrostatic difference 
between Sr 2+ and La 3+ . The long-range Coulomb inter- 
action is characterized by the dielectric constant en, he., 



the Coulomb energy Eq 



/er>a. The model is de- 



scribed in terms of a generalized Hubbard Hamiltonian 



H = Hv 



H , 



(1) 



with the ordinary Hubbard Hamiltonian including 
nearest-neighbor hopping and the on-site repulsion 



n iT n a 



and a part taking account of the long-range electron-ion 
and electron-electron Coulomb interaction 



= E^+5£ 



h z Wijhj, 



with the interaction matrix = Ec/\ri — Tj\. For 
charge neutral systems, the model considered here 
[Eq. involves three different parameters, namely 

Ec/t, N and U/t. The relative strength of the long- 
range Coulomb interaction Ec/t mainly affects the self- 
consistent screening of the ions. A realistic value for 
the experimental system realized by Othomo et aZ.— is 
Ec = 0.8t corresponding to £d ~ 15, a ps 3.9 A and 
t w 0.3 eV^^ For the following discussion we fix Eq 
at this value. The relative on-site repulsion U jt is in- 
trinsically large for Mott insulators. Nevertheless we will 
treat this parameter continuously and consider the case 
where magnetic order is suppressed - therefore tuning the 
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central material from a paramagnetic metal to a param- 
agnetic Mott insulator. Furthermore, we focus in our 
work on a systematic study on how the T — properties 
change by varying the width N of the Mott-insulating 
material. 



B. Mean-field treatment 

One way to arrive at an effective low-energy theory 
is provided by the KR SBMF approximation. This can 
be done in full analogy to the homogeneous single-band 
Hubbard model. 17 We thus reformulate the model Hamil- 
tonian [Eq. ([I])] in an enlarged Hilbert space by introduc- 
ing pseudo fermionic operators /£P and auxiliary bosonic 

operators (slave bosons) ej , p^J and standing for 
the creation (annihilation) of empty, singly occupied sites 
with spin a and doubly occupied sites, respectively. Yet, 
this new Hilbert space is much larger than the Hilbert 
space of dJ). To project on the physical subspace we need 
to fulfill the following constraints for each individual site 
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$>U,+44 = J. (2) 

(T 

Pljiv + fidi = flL- (3) 



In the mean-field approximation the slave-boson fields 
are replaced by their classical values and the constraints 
[Eq. (|2|3| ] arc fulfilled on the average. We will concen- 
trate in the following on the paramagnetic mean-field so- 
lution. This allows to discuss important aspects of the 
Mott physics inherent in the model. In what follows we 
will use the convention i = (j, I) where I is the layer in- 
dex (z-direction) and j specifies the site in each individ- 
ual layer (x, y-direction). We fix the origin z = in the 
center of the Mott-insulating material and set the lattice 
spacing a = 1. As independent mean fields we choose 
the amplitude of doubly-occupied sites di and the total 
charge distribution ni in layer I as well Lagrange 
multiplier A; which controls the self-consistency of the 
charge distribution. The mean-field Hamiltonian can be 
diagonalized by a canonical transformation 



filer 
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which is composed of an in-plane Fourier transformation 
and of a transformation to the solutions of an effective 
one-dimensional Schrodinger equation in the z-direction 



(z, 2 -l)e fc + Aj <M0 



t 2J ZlZl+^lpkvil + 7) 
7=±1 

{E vk -£ k )il> k „(l). (4) 



The quantum number v = 1, 2, . . . is the band index, i. c. 
it labels the solutions of Eq. (0| for fixed k. The bound 
states of Eq. ^ define bands of a quasi-two-dimensional 



metal. Furthermore, the symmetry z <-> — z allows to 
choose eigenstates with a fixed parity (— l)" -1 , 



lM-0 = (-1)"-Vkv(f). 



(5) 



The in-plane hopping is renormalized by a factor zf - 
corresponding to a layer dependent mass-rcnormalization 
ml jm = 1/zf - and the hopping between different layers 
by zizi' wherein 



y/(l -m+ df)(ni - 2d'f) + diy/ni - 2d'f 
zi(ni,ai) = == . 

(6) 

In Eq. ([!]), the inhomogeneous setup leads to an effec- 
tive one-dimensional potential (zf — l)s k + A; along the 
heterostructure which can vary essentially for different 
values of the transverse momentum fc, i.e., the in-plane 
kinetic energy e = £k = — 2t(cosk x + cosk y ). The diag- 
onalization of Eq. ((4|) therefore has to be performed for 
each Ek separately. 

The mean fields di, n; and A; are determined by the 
stationary point of the free energy density at the inverse 
temperature {3 = 1/T, 



f(d,n,X) 



) 



(7) 



+U T Y, $ + \ E mWwn v - ^(A, - Vr)n h 
1 u> 1 

under the constraint of charge neutrality ni = 
The saddle-point equations incorporating this condition 
arc nonlinear, coupled self-consistency equations for the 
mean-fields: 



N. 



m - ni, 

rr J dz l 

odi 



t E 2 H 



ddi 



(8) 
(9) 



7=±1 

dz 2 



9/ V 7, ^ V (7) 



7=±1 

where we have introduced the following definitions 
2 
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I kv 



(10) 

(11) 

(12) 



ku 



Xi 



(7) 



— Y,^u{l)M + l)fT{E kv ), (13) 



with the Fermi function f T {w) = (1 + e^) -1 . The 
saddle-point equations (|5j)- (fTTJ]) together with the effec- 
tive Schrodinger equation ^ are the basic equations in 
our approach. 
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We remark here that in the physical subspace the op- 
erator identity 

Uh^hii = (U - U )d\d i + Uoh^hii (14) 

holds with an arbitrary constant Uq. However, the mean- 
field results depend on the choice of Uq. The freedom 
resulting from Eq. (|14[) is not unique to the particu- 
lar model considered here but for example also exists 
in the SBMF approximation of the homogeneous Hub- 
bard model. However, in the latter case we demand that 
for [7 = the non-interacting case is reproduced by the 
SBMF approximation^ and consequently Uq = 0. In the 
model considered here we find it useful to choose a finite 
threshold Uq = Eq. Because usually the physical rele- 
vant parameter regime is U > Eq we define a minimal in- 
teracting limit by U = Eq and demand that in this limit 
the SBMF approximation reproduces the Hartree mean 
field approximation. Hence, we set Uq = Eq meaning 
that double occupancy is only significantly reduced for 
U > Eq. In our approach the relevant parameter describ- 
ing the correlation induced deviation from the Hartree 
mean- field results is therefore U T = U — Eq ■ 

Including the shift for the on-site term discussed above 
the long-range Coulomb interaction matrix in units of Eq 
is given by 

Ww_ = (l-Sw) 

E c Jn 2 + m 2 + (l- I') 2 

I 7+ E / / , I • ( 15 ) 

Since we consider only charge neutral systems the diver- 
gent part of Eq. (jTSJ) is canceled by the divergent part of 
the long-range ion-electron interaction. 
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FIG. 2: (Color online) The mean-field d, for the N = 10 ion- 
layers heterostructure for Ur = 0, . . . , 34t. The fraction of 
doubly occupied sites in layer I is given by df . 

Here, K(x) is the complete elliptic integral of the first 
kind. Note that there is a van-Hoove singularity at e = 0. 
For the diagonalization of Eq. ^ we impose open bound- 
ary conditions and take typically about 30 layers into ac- 
count, depending on the number of ion layers N. All the 
calculations are performed at T = 0. The fact that only 
the band energy e enters formally into the SBMF analy- 
sis means also that we can treat various lattice topologies 
in the exactly same framework with modified density of 
states. 



C. Numerical scheme 

One way to solve the nonlinear saddle-point equations 
is first maximizing the free energy f(n, d, A) with respect 
to A and then minimizing the resulting function with re- 
spect to d and n (see, e.g., Ref. H^). It turns out that 
this strategy leads to a stable algorithm in a wide range 
of parameters whereas an iterative solving of the saddle- 
point equations fails for larger values of U due to the 
strong nonlincarity. Thus, our algorithm contains basi- 
cally two loops where in the first loop the free energy 
is maximized with respect to A and in the second loop 
the resulting function is minimized with respect to d and 
n. The calculation of the free energy and of its gradi- 
ents in the thermodynamic limit Ni\ — > oo involves a 
two-dimensional fc-intcgration which can be reduced to a 
one-dimensional energy-integration including the density 
of states of the two-dimensional nearest neighbor tight- 
binding model 

N(s = f (l/2n 2 t)K (1 - ( £ /4i) 2 ) if |e| < 4*; m 
' I else. ^ ' 



III. GROUND-STATE PROPERTIES 
A. Double occupancy and charge distribution 

We first discuss the fraction of doubly occupied sites 
given by the mean field df and the charge distribu- 
tion m. It is clear, that in order to reduce the on- 
site potential energy cost, doubly occupied sites are sup- 
pressed by increasing f/ r . We show this in Fig. [2] for the 
N = 10 heterostructure. In the minimal interacting case, 
U t = 0, the fraction of doubly occupied sites in layer I 
is df = nf/4. With increasing U T , di tends to zero. As 
expected, doubly occupied sites are strongest suppressed 
in the Mott-insulating material, namely in the center of 
the heterostructure where the charge is close to 1 (sec 
Fig. [3]). Consequently, the hopping renormalization fac- 
tor z\ [Eq. ([6])] is reduced there reflecting the suppression 
of charge fluctuations. On the other hand, in the band- 
insulating material, di varies only little with increasing U r 
and stays close to the minimal interacting value. There- 
fore, zi w 1 in this region. 

In the paramagnetic phase discussed here, the system 
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FIG. 3: (Color online) Charge distribution for the TV = 10 
ion-layers heterostructure for various values of U? between 
and 34t. 



has two competing possibilities to reduce the double oc- 
cupancy and thus the on-site energy cost: 

(i) Spreading of the electron charge distribution. 

(ii) Localization of the electronic states. 

Of course, in both cases, long-range Coulomb as well as 
kinetic energy has to be payed. Figure [3] shows the elec- 
tron charge distribution rii of the TV = 10 heterostructure 
for different values of U r . We find a rather small depen- 
dence of ni on U r : the effect of the long-range Coulomb 
field dominates over the short-range correlations and con- 
sequently the above mentioned two mechanisms play a 
rather minor role for the overall growth direction electron 
density distribution.— ' 14 ' 15 Nevertheless, a closer look at 
the evolution of the charge distribution allows to distin- 
guish two different regimes characterized by one of the 
two mechanisms. Starting from the minimal interacting 
case and increasing U T results in a spread of the parti- 
cle density. This effect is also included in the Hartree 
approximation where the decoupling of the U -term leads 
to an additional term in the self-consistent electrostatic 
potential A; — > A; + U T ni/2 which favors a homogeneous 
charge distribution. We therefore refer to this regime as 
the Hartree regime. But for all TV > 1 we observe that 
in most layers the density ni evolves non-monotonically 
with increasing U r . Namely, there is a turning point 
where the slope of ni(U r ) changes sign. Figure |4] shows 
the occupancy of the central layers I = 0, I = ±1 and 
I = ±2 of the TV — 10 heterostructure as a function of 
U T . The turning-point- value of U r is different for the dif- 
ferent layers and varies also with TV but is roughly in the 
range of the critical value L^? R ~ 16t of the Brinkmann- 
Rice transition in the homogeneous system. We therefor 
refer to the regime U r > Uf K as the Mott regime. 

These observations suggest to look at the following two 
quantities characterizing the charge distribution. We de- 
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FIG. 4: (Color online) The occupancy of the central layers as 
a function of £/ r for TV = 10. 
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FIG. 5: The mean distance P of the charge from the center 
as defined in Eq. (|17p shown for TV = 3, 5, 8, 10 as a function 
of U T . 



fine the square mean distance of the charge from the cen- 
ter of the heterostructure as 



N 



(17) 



and a quantity measuring the width of the interface as 
A = 53nj(l-n,). (18) 

A is a measure of the region where the electron density 
ni differs significantly from the bulk values n; = (band 
insulator) and ni = 1 (Mott insulator). In the strongly 
interacting limit the mobile carriers are mainly confined 
to this region (see also Sec. IIII E|) . Figure [5] shows I 2 
for TV = 3, 5, 8, 10 as a function of U r . In all cases we 



FIG. 6: (Color online) The quantity A = ^ n;(l — ni) mea- 
suring the width of the interface as function of U r for different 
heterostructures N. 

can clearly distinguish Hartree and Mott regime. In the 
Hartree regime we observe a rather steep increase of I 2 
reflecting that the reduction of the doubly occupied sites 
is mainly achieved by the spreading of the charge, namely 
by the mechanism (i). On the other hand, in the Mot 
regime, I 2 is approximatively constant (it flattens or even 
slightly decreases) indicating that in this regime the effect 
of localization is predominant [mechanism (ii)]. A similar 
behavior as function of U T is also observed in the width 
A of the interface (see Fig. [5]). In addition, for N > 5, 
A becomes almost independent of N and differs only in 
the vicinity of R reflecting that for large N the two 
interfaces are disconnected. 

The behavior of I 2 and A is a consequence of the com- 
petition between the three energy scales in the problem, 
namely Ec, U r and the bandwidth W = 12t. In the 
Hartree regime Eq < U T <C W the qualitative behavior 
observed in Fig. [5] and [5] is basically determined by the 
competition between Ec and C/ r . The penetration depth 
thus increases by increasing U T due to the mechanism (i). 
On the other hand, in the Mott regime Ec < W <C U T 
the penetration depth is mainly determined by the com- 
petition between Ec and W and therefore essentially in- 
dependent of U Y . Because the electronic states are almost 
localized or, in other words, the binding of the electrons 
to the ions is strengthened, the charge distribution is al- 
most independent of U r . 

B. Renormalization of the quasiparticle dispersion 

We now go over to a discussion of the solutions of the 
effective Schrodingcr equation ^ . Whenever it is appro- 
priate, we will use the notation E v (e) for the quasipar- 
ticle dispersion and ip ev (V) for the transverse part of the 
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FIG. 7: The quasiparticle dispersion Eku of the partially filled 
bands in the N = 3 and N = 8 heterostructure as a function 
of Ek = — 2t(cos k x + cos k y ) both for U T = and U r = 24t. 

quasiparticle wave function because the momentum de- 
pendence enters only via e = Ek = — 2f(cos fc^-l-cos k y ). In 
general we find N + 5 partially filled bands in the ground 
state. In the minimal interacting limit (U T = 0), the dis- 
persion of these bands has the form E v (e) = E v + e, i.e. 
E v (e) has slope one. Increasing U T leads to a renormaliza- 
tion of the quasiparticle dispersions and E v {e) flattens, 
thus indicating localization of the electronic states. The 
slope of E v {ek) is given by 

Z kv = ^=Y^zUUlf (19) 

which, as a consequence of the layer-dependent mass 
renormalization, is in general less than one. 

As an example we consider here the N = 3 heterostruc- 
ture with 8 partially filled bands and the N = 8 het- 
erostructure with 13 partially filled bands. The disper- 
sion in the minimal interacting case and for U r = 24t is 
plotted in Fig. \7\ The particular geometry of the sys- 
tem, namely the symmetry z <-> —z, allows that some 
of the quasiparticle bands are (almost) degenerate. We 
can distinguish two types of degeneracies: (1) Essen- 
tially independent of the interaction strength we find that 
the almost unbounded states (the highest partially filled 
bands) give rise to pairwise degenerate bands. In the case 
N = 3 this is observed for the states v = 5, 6 and v = 7, 8. 
A similar feature can be seen for the N = 8 heterostruc- 
ture where the six highest bands v = 8 — 13 are pairwise 
degenerate. (2) In addition, by increasing J7 r , also the 
lowest bands can become pairwise degenerate in certain 
fe-regions. This is observed in Fig. [7] for the ^ = 1,2 and 
v = 3, 4 bands of the N — 8 heterostructure and reflects 
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FIG. 8: (Color online) The quasiparticle dispersion near 
the Fermi energy for = 2At. Two types of avoided level 
crossing are observed depending on whether N is even or odd. 
The type of avoided level crossing can be understood within 
a minimal model considering 3 (N even) respectively 4 layers 
(TV odd). 



that the two interfaces become disconnected from each 
other. The aspect of correlation induced degeneracy will 
be discussed in details in the next section. 

The renormalization of the quasiparticle dispersion is, 
in general, the stronger the bigger N and C/ r and can 
substantially depend on v and e = A qualitative 
understanding of the fc-dependence is given by the fol- 
lowing considerations. For the moment let us ignore the 
interlayer hopping and assume that the bulk system of 
the sandwiched material is Mott-insulating. Then the 
dispersion in each individual layer can be interpreted as 
an uniformly rcnormalizcd band e;(e) = zfe + A;. Al- 
most entirely flat bands near the Fermi energy, which 
are close to half filling, correspond to layers in the center 
of the heterostructure. On the other hand, more disper- 
sive bands with low filling correspond to layers further 
away from the center. If we now allow for a finite inter- 
layer hopping Viv = —tz\Z\> the uniformly renormalizcd 
bands hybridize. As a result, the renormalization of the 
hybridized bands can strongly depend on e = 

The fc-dependency is strongest in the vicinity of the 
Fermi energy where E v (e) shows a complex pattern of 
avoided level crossings as function of e (see Fig. |5J). De- 
pending on whether N is even or odd we can distinguish 
two types of avoided level crossing. The "building blocks" 
can be calculated within a minimal model including 3 (N 
even) or 4 (N odd) layers as shown in the Appendix IA1 

The above discussion suggests a certain relation to 
the heavy-fermion ground state of the Periodic Anderson 
models There, the hybridization between the disper- 
sionless and strongly correlated / level and the itinerant 
conductance band leads to a heavy Fermion liquid with 
strongly rcnormalizcd quasiparticle properties 3i How- 
ever, similar to the periodic Anderson model, in the 
model considered here the stability of the paramagnetic 
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FIG. 9: The effective one-dimensional potential A; + {zf — l)ek 



for various values of 
N = 10. 



-4i, -3i, 



, it for U r = 14t and 



heavy-Fermi liquid state against magnetic ordering is an 
important question which will be addressed elsewhere. 

Finally we remark here that the filling of the quasi- 
particle bands also changes by varying U r . However, it 
seems that this effect is rather weak and we refrain from 
discussing this point in more details. 



C. Correlation induced degeneracy and momentum 
dependence of tpkv(l) 2 

Within the SBMF approximation the fc-dependence 
of the quasiparticles formally enters through the term 
(zf — l)ek in the effective Schrodingcr equation This 
leads to a e-dependent effective potential A; + (zf — l)e 
which includes the screened ion-potential as well as a 
'Mott potential' introducing an effective central barrier 
(— » double-well potential form) for e < and a single 
well for e > 0. Thus, for fc-values in the center of the 
first Brillouin zone (e < 0), a double- well potential is 
formed and the regions I > and I < arc almost decou- 
pled. As an example we show A; + (zf — l)e in Fig. [9] for 
N = 10 and U r = 14i. This means that electrons which 
have a low (in-plane) kinetic energy prefer to stay mainly 
in the interface region rather than in the center where the 
effect of the on-site interaction U r is much stronger. To il- 
lustrate the rather complex dependence of the transverse 
part i/j eu on e, we show in Fig.[TO]a density plot of ipev(l) 2 
for v = 1, . . . , 4 as a function of e and I for the same pa- 
rameters as in Fig.[9j In this plot we interpolate ^^(i) 2 
between different Vs. The e value of the Fermi energy 
is denoted by e* and is indicated as a white dotted line. 
Due to the double-well structure the density plot of tpei 
and looks very similar for e < 0. As a consequence, 
energetically nearby states of different parity can become 
almost degenerated in a certain e-range. Figure [TT] shows 
the four energetically lowest lying bands of the = 10 
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FIG. 12: Quasiparticle weight Z v of the partially filled bands 
v = 1, .... 8 as a function of U T for TV = 3. 



FIG. 10: (Color online) Density plot of the square of the four 
lowest lying transverse wave functions Tp elJ (l) 2 as a function 
of I and e for the same parameters as in Fig. [5] The Fermi 
energy is indicated as white dashed lines. 




(7T, 71") 



FIG. 11: The quasiparticle dispersion Eku of the four lowest 
lying bands in the N — 10 heterostructure for U T = 14t. Near 
(k x ,k y ) w (0,0) they are pairwise degenerated. 



heterostructure for U r = 14t. These bands arc pairwise 
degenerated near (k x ,k y ) ~ (0,0). 



D. Quasiparticle weight 

A characteristic quantity describing correlation effects 
is the quasiparticle weight Z v of the band v. It for exam- 
ples enters the linear specific heat at low temperatures 



c v N 



-T 



Q T 2 



7T 

~3~ 



Z u 



T, 



where E v (e*) = and N(e) is the free electron density of 
states [Eq. (fl"6]) ]. The quasiparticle weight accounts for 
the renormalization of the Fermi velocity - a vanishing 
Z v thus indicates a localized state. It is given by the 
value of Zku [Eq. (|19| ] at the Fermi energy: 



Z„ 



8Ey 

de 



(20) 



Here, ip e * v is the transverse part of the wave function 
at the Fermi energy. The above expression will reappear 
later when we calculate the discontinuity in the momen- 
tum distribution function at the Fermi energy at T = 0. 

The quasiparticle weight Z v of the partially filled 
bands of the N = 3 heterostructure is shown in Fig. [12] 
as a function of J7 r . Energetically lower lying bands are 
in general more affected by U r and the v = 1 band has 
always the smallest quasiparticle weight. Note however 
that this band is not closest to half filling. The exact or- 
der of Z v for v > 1 depends on N and in some cases also 
on U T . Because odd wave functions have nodes at z = a 
parity effect is visible. If N is even the electronic sites are 
at I = 0, ±1, ±2, . . . and odd-parity bands are less renor- 
malized because ijjetv{0) = 0- On the other hand, if N 
is odd, the electronic sites are at I = ±1/2, ±3/2 . . . and 
the odd-parity bands tend to be more strongly renormal- 
ized because the zero at z = lies between the electronic 
sites. For N = 3, this parity effect is seen in Fig. [T2l where 
Z4 < Z3 and Zq < Z5. We note here that both mag- 
nitude and order of the quasiparticle weight is in good 
qualitative agreement with the DMFT results shown in 
Fig. 3.(b) of Ref.0 

As a further example we show in Fig. [13] the quasi- 
particle weight Z\ of the lowest lying band as a func- 
tion of U r for various values of N. Note that for = 1 
there is no significant renormalization of the quasiparticle 
weight which is in agreement with the observation made 
in density-functional theory based calculations for the 
N = 1 heterostructureii and reflects the special role of 
this case 7 . In order to have a more quantitative descrip- 
tion how on-site correlation effects increase with growing 
N we define a characteristic value U^N) for each N by 
extrapolating Zi(U r ) at the inflection point to zero which 
yields U^N). In Fig.[Tl]we show U^(N) as a function of 
N. For JV-^oowe find that C/ r c saturates at [/ C BR w 16i, 
the critical interaction strength for the Mott transition 
in the bulk system of the sandwiched material. 

Beside the suppression of z; in the center of the het- 
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FIG. 13: The quasiparticle weight Z\ for N = 
1, 2, . . . , 8, 10, 12, 15, oo as a function of U T . In order to de- 
fine a characteristic value U^(N) we extrapolate Zi at the 
inflection point to zero. 
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FIG. 15: (Color online) The transverse part of the wave func- 
tion i\>siv(V) at the Fermi energy for the two lowest lying bands 
for various values of U r = — 34t for the N = 10 heterostruc- 
turc. 



on-site energy cost. Among all the states characterized 
by (kv) only a part is strongly affected by the correla- 
tions. This allows for an optimal reduction of the on-site 
energy cost without cutting back much from the kinetic 
or long range Coulomb energy. Thus, the possibility that 
correlation effects can be selectively split among a large 
number of different states prevents the system from un- 
dergoing a strict Mott transition. 



FIG. 14: I!? as a function of N. The dashed line denotes the 
value of the Mott transition in the TV — > oo model. 

erostructure there is also a slight change in ip £ * u with 
increasing [/ r . Namely, there is a clear tendency towards 
localization in the center of the heterostructure which re- 
duces the quasiparticle weight Z v [Eq. (PO)) ] even more. 
As an illustrative example, the wave functions at the 
Fermi energy of the N = 10 heterostructure are shown 
in Fig. [15] for v = 1, 2 and for a range of U r = 0, . . . , 34t. 

Finally, we point out that we do not observe a strict 
Mott transition for finite N and U T in the sense of a van- 
ishing quasiparticle weight Z v respectively of a divergent 
effective mass m\ ~ Instead, the mutual doping 

of Mott and band insulator leads to metallic behavior of 
the whole structure. The system gains energy by keep- 
ing a certain amount of charge fluctuation meaning that 
Z v > 0. In other words, the hopping between the indi- 
vidual layers is never renormalizcd to and n\ < 1 and 
z\ > 0. In the spatially nonuniform system discussed 
here, there is a large amount of freedom to reduce the 



E. Layer resolved spectral density and coherent 
particle density 

The dynamical properties of the correlated electrons 
in the heterostructure can be best studied by looking 
at the single particle spectral function. In principle, 
the spectral functions are accessible in photocmission^S 
or scanning tunneling microscopy. In the framework of 
the SBMF approximation only the coherent (low-energy) 
part of the spectral weight can be captured. In or- 
der to access to lower and upper Hubbard bands within 
slave-boson theory one has to go beyond the saddle-point 
approximation^ which is beyond the scope of this pa- 
per. Nonetheless, the low-energy part is important for 
the metallicity of the heterostructure. 

The fc-depcndcnt layer resolved spectral density is de- 
fined as the imaginary part of a retarded one-particle 
Greens function 

4<7(w)=--ImG (to (w). (21) 



10 



Within the SBMF approach the coherent part of Guka (u>) 
is given by 



Gcoh 
Ilka 



(«) = £ 



(22) 



E kv + 10+ 

and we find for the coherent part of the spectral density 
O) = £ aftM0 2 % - (23) 

Integration of A;fco-(w)/o(a;) over l,ui yields the ground- 
state momentum distribution function of the correlated 
electrons. The coherent part of the momentum distribu- 
tion function n£° h = contains contributions of 
the different bands 



^coh ^coh 



b kva 



L kv 



Zkufo(Ekv) 



(24) 



where Z kl , is given in Eq. Q19p. This means that each 
state (kva) is weighted by a factor Z kv < 1. There- 
fore, also the discontinuity of at the Fermi energy 
is reduced, namely, it is equal to Z v = zfip E * v {l) 2 = 
dE u (e*)/de and we recover the expression for the quasi- 
particle weight given in Eq. (f20l) . 

The (k,a integrated) layer resolved spectral density 
can be expressed in terms of the free electron density of 
states N(e) [Eq. I|16p] and is given by 



1 1 k<r 



(25) 



= lo and 



Here, e ULU is defined by the equation E u {e u 
Z yuJ = dE v (e VUi )/de. 

Figure [TBI and [T71 show the layer resolved spectral den- 
sity A™ h (uj) in the minimal interacting case and for 
U r = 24i, respectively, for the N = 10 heterostructurc. 
In the latter case, the coherent part is strongly reduced 
towards the center of the heterostructure. The missing 
spectral weight indicates nearly insulating behavior and 
that most of the spectral weight is transferred to up- 
per and lower Hubbard band at a higher energy scale. 
(Note that the SBMF approximation is not able to catch 
the incoherent spectral weight which defines upper and 
lower Hubbard band.) However, outside the heterostruc- 
ture, the spectral weight of the near-Fermi-surface-states 
is only weakly suppressed and almost identical to the 
minimal interacting case. We find that in the band- 
insulating material the chemical potential (u> = 0) is, in 
all cases studied, pinned to the the bottom of the conduc- 
tion band. The highest-lying electron states are therefore 
only weakly bound (see also Ref. 0). The various sin- 
gularities showing up in the spectral density enter here 
because of the van-Hoove singularity in the free electron 
density of states iV(e) in combination with the infinite 
life-time of the quasiparticles at T = 0. At finite temper- 
ature we expect that these singularities are washed out 
due to finite-life-time effects ~ T 2 of the quasiparticles. 
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FIG. 16: Layer resolved spectral density for U T = for the 
N = 10 heterostructure showed for Z = 0, .... 9. 
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FIG. 17: Layer resolved spectral density for U T 
N = 10 heterostructure showed for Z = 0, .... 9. 
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2U for the 



By integration of Eq. ([25)1 we find the coherent part of 
the particle density 



^coh 



dw^ oh (w)/ T M = zfm 



(26) 



Figure [18] shows 



-.coh 



for various values of the on-site 



interaction U T . It is interesting to note that the ratio 
nf oh /ni = zf = m* /m gives the mass renormalization in 
layer I. A similar observation was mentioned in Ref. 
on the basis of quasiparticle dispersion calculations in 
the framework of DMFT. The coherent particle density 
clearly shows humps at the interfaces between Mott and 
band-insulating material whereas in the center of the het- 
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FIG. 18: (Color online) The coherent particle density n^ oh 
of the TV = 10 heterostructure for various values of U T = 
OA, . ..,34t. 



erostructure the coherent part is strongly suppressed for 
strong on-site interactions. These humps have a width 
of approximately 3a which suggests that the coherent 
quasiparticles responsible for the metallic behavior are 
confined to a very narrow region at the interface. For 
large TV and U r the coherent particle density is essen- 
tially independent of TV and U T . 



F. Drude weight and optical conductivity 

The zero frequency contribution to the optical conduc- 
tivity at T = - the so called Drude weight - is one of the 
most relevant quantities both experimentally and theo- 
retically to distinguish between an insulator and a metal. 
It was pointed out by Kohn^i that for a system with 
a (discrete) translation invariance in e.g. n-direction the 
real part of the complex optical conductivity a = a' + 
has the forrn^ 



a'(u) = D5{uj) 



(27) 



with lim^^o CJ<7 reg( CJ ) = 0- The weight of the Dirac delta 
function S(cu) is called the Drude weight D. D provides 
an order parameter for the Mott transition at T = 
in non-disordered systems^^ The system is metallic if 
D > and insulating if D = 0. Studying the Drude 
weight therefore allows for an alternative assessment of 
the metallicity of the interface between band- and Mott- 
insulating materials. 

As we have shown in detail in the previous sections the 
low-energy behavior of the system is described, at least 
within the SBMF approximation, by coherent quasipar- 
ticles (see also the discussion in the last section of this 
paper, Sec. IIV|) . In this case, also the low-energy part of 
<t{ui) is expected to be dominated by the contributions of 
these quasiparticles^ A reasonable way to estimate the 



low-energy behavior of the optical conductivity is there- 
fore to apply linear response theory to the renormalized 
mean-field Hamiltonian. In this way a close formal anal- 
ogy to previous studies of the optical conductivity in the 
single-band Hubbard model can be &v&wnM£2£^&2££L 
Alternatively, the Drude weight D can be obtained from 
the ground-state energy of a system with twisted bound- 
ary conditions - an idea originally formulated by Kohn. 31 
In the Appendix [B] we show that the same result for the 
Drude weight is directly obtained by considering twisted 
boundary conditions and using the SBMF approximation 
in order to estimate the dependence of the ground-state 
energy on the twist. 

In the following we choose units such that h = 1 = c. 
The response of the system to a time dependent uniform 
electric field in the x-dircction is found to be 



uj + i,0+ 



(28) 



The kinetic energy of the coherent quasiparticles 



kv 



is given by the expectation value of the renormalized ki- 
netic energy operator in ^-direction 



Tr 



-2t^2 cos k x zffl a f lko 



klo 



The current-current correlation function is given by 



H = ]v"^l <ri|Jp:c|0> 



i 



u - (E n - E ) + iOr 
1 



to + [E n - E ) + i0+ 



(29) 



with the renormalized paramagnetic current in in- 
direction 

J P x = e2t sin k x zf f lka f lkcr . 

klcr 

The many-body states {\n}} in Eq. (f29|) arc pscudo- 
fermionic Slater-determinants with energy E n . We de- 
compose x m t° real and imaginary part, x = x' + *x"- 
The real part of <j(uj) is found to be of the form of Eq. (|2"T|) 
with the Drude weight given by 



D 



tie 



-T 



7TX'(0) 



(30) 



and a regular part a' zeg (u>) = —x" { UJ )/ LU - The first term 
in Eq. (|30p is positive and converges quite rapidly both 
with increasing U r and N because for large N and U T the 
coherent quasiparticles are located in a narrow region at 
the interfaces between the two materials, as discussed in 
the previous section. The rapid convergence is shown 
in Fig. [T9l where the TV-dependence of (— T x )q for fixed 
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FIG. 19: The Drude weight [Eq. ((30]) ] and its two quasiparticle 
contributions as a function of TV for U T = 26t. 



FIG. 20: Low-energy optical conductivity of the TV = 3 
haterostructure in the strongly correlated case Ur = 24i. The 
main contributions to the different features can be assigned 
to transitions between two bands v <-> v' and are shown as 
thin lines. don't forget to add the prime... 



U T = 26i is plotted. The second term ttx'(O) in Eq. JSOj) 
is negative, thus it reduces the Drude weight. However, 
this weight is transferred to the regular part <7 reg at finite X (0) 
values of co. This is essentially a consequence of the /- 
sum rule for the real part of <t(w) 



a'(uj)du; = ——(-f x } 



2Nn 



(31) 



meaning that 



a ' {u>)du). 



In general x{lo) ^ because it contains non-vanishing 
contributions from matrix elements between band pairs 
of the same parity 



{vk\J px \v'k) = e2tsmk x } j z l i/)kv(l)i>kv'(l)- 

al 

Hence, the fact that the renormalization factor z; de- 
pends on I formally allows for interband transitions 
when a time-dependent electric field is applied in the in- 
direction. However, this is only possible for U T > - 
in contrast to the case of an electric field applied along 
the z-dircction.— As shown in more details in the Ap- 
pendix[B] this is closely related to the fact that the renor- 
malization of the quasiparticle dispersion Z kv [Eq. (fTT))) ] 
is fc-dependent. Therefore, unlike the case of the uni- 
form Hubbard model, in systems with spatial inhomo- 
geneity x(w) acquires non- vanishing contributions from 
interband quasiparticle scattering processes when U r > 0. 



For the real part at u> = we explicitly find that 



2 v \(n\J px \0)\ 2 
iV,,^ E n -E 



\(vk\J px \ V 'k) 



2 fo(E kv >) - fo(E kv ) 



kv' <iv 

^E^i^i 2 /o(^). 

I kv 



Ekv — E>kv> 



(32) 



We expect that also x'(0) saturates with increasing U T 
and N because in the large TV and strongly interacting 
limit essentially two independent interfaces contribute to 
the Drude weight D. However, as seen in Fig. Q1J1 the 
convergence of x'(0) with increasing for fixed U r is 
much slower than for (—T x )q. The regular part of the 
low-energy optical conductivity a' Teg (uj) = — x"(w)/ui can 
be written as 



2tt 

fo(E kv )]6(u; 



J2 \{vk\j px \ v 'k) 



I kv'<u 
2 



{E kv - E kv ,f). (33) 



Figure [2lH shows <r' Icg (uj) of the TV = 3 heterostructure in 
the strongly interacting case U T = 24t as calculated from 
Eq. (|33j) by approximating the (5-function by a Lorentzian 
of width O.Olt. The main contributions to the different 
features can be assigned to transitions involving only two 
different bands u <-> v' of the same parity. The width of 
the individual features is due to the fc-dependent splitting 
of the bands and not due to the approximation of the 5- 
function. 

The £/ r -dependence of D is shown in Fig. [21] for differ- 
ent heterostructures of width TV. Again we can distin- 
guish between the Hartree regime and the Mott regime 
that are clearly separated by a common intersection point 
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FIG. 21: Drude weight for different values of AT as function 
of U T . 
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FIG. 22: The maximal value of the Drude weight and the 
corresponding N. 



of D(N, U T ) slightly below U T « C/ C BR for AT > 5. This 
two regimes are characterized by the different behavior 
of D in the limit N — > oo: In the Hartree regime, D in- 
creases with increasing N because the number of mobile 
carriers is proportional to N. In the Mott regime the AT 
dependence is much weaker and the main contribution 
to the Drude weight stems from the interface regions. D 
decreases with growing N and we expect that it saturates 
eventually. 

Because the dependence of D on U T is rather weak if 
N is small (1 ~ 4) there is another interesting aspect. 
Namely, in the Mott regime there is a finite value of N 
which maximizes the Drude weight for a given interaction 
strength. We show in Fig.[21]-Dmax and N max as function 
of U t . Both are decreasing functions of U T . 



IV. CONCLUSIONS 

To summarize, we have presented a mean-field study 
of a model heterostructure characterized by a BI/M/BI 



stacking (see Fig. [T]) based on the Kotliar-Ruckenstein 
slave-boson theory. Although the model discussed here 
is too simple to make quantitative statements about the 
experimentally realized LaTiOa/SrTiOs interface it al- 
lows to discuss some aspects of short-range correlation 
in spatially nonuniform systems. In particular, it sug- 
gests that the mobile carriers are confined to a relatively 
narrow region at the interface^ forming a quasi-two- 
dimensional electron gas as experimentally observed in 
LaTi0 3 /SrTi0 3 heterostructures£ 

The SBMF approximation offers a relative simple and 
transparent way to include both the self-consistent deter- 
mination of the charge distribution as well as spatially 
varying hopping renormalization factors due to strong 
correlations. The basis of this approach is an effec- 
tive one-dimensional Schrodinger equation in combina- 
tion with self-consistency equations for the charge dis- 
tribution and the double occupancy. Due to the lack of 
spatial homogeneity the quasiparticles obtain a nontriv- 
ial momentum dependence which can vary substantially 
by changing the on-site interaction U or the width N 
of the sandwiched material. This dependence was stud- 
ied by calculating the renormalization of the quasiparti- 
cle dispersion and the quasiparticle weight. Dynamical 
properties were studied by calculating the layer resolved 
spectral density and the low-energy part of the optical 
conductivity including the Drude weight. To which ex- 
tend the quasiparticle properties discussed in this paper 
are reliable depends on the interesting question whether 
or not the quasi-two-dimensional electron gas formed at 
the interfaces is a Fermi liquid (see also the discussion in 
Ref. H). Whereas Sri-^La^TiOs shows clear Fermi liquid 
behavior with enhanced electron-electron scattering pro- 
cesses as x — ► 1 (Ref. [H) and resistivity measurements 
on LaTiOa/SrTiOs supcrlattices are consistent with this 
picture^ it is an open problem if the Fermi liquid de- 
scription survives when reducing the dimensionality fur- 
ther to a quasi-two-dimensional system. However, in the 
SBMF analysis of the model, the system is a Fermi liq- 
uid with well defined quasiparticles and the present ap- 
proach offers a rout to access to Fermi liquid parame- 
ters on the basis of a microscopic model. On the other 
hand, we confirmed various results reported in previous 
dynamical-mean-field studies^ which leads to the con- 
clusion that the SBMF approximation gives a reasonable 
description of the low-energy part and is thus able to 
catch the most important aspects of correlations in spa- 
tially nonuniform situations. However, only on-site cor- 
relations can be taken into account and the introduction 
of a layer dependent mass renormalization distinctly in- 
creases the computational effort as compared to Hartree 
mean-field calculations. But for all that, compared to 
DMFT calculations, the SBMF approximation has the 
advantage that it is numerically less involved. 

So far we did not treat other important aspects inher- 
ent in strongly correlated electron systems. The most 
conspicuous is the ignorance of magnetic ordering which 
especially in bipartite lattices is of relevance. To some 
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extent this can be included also on the mean-field level 
by considering an enlarged unit cell. On the other hand, 
the results shown in this paper are basically independent 
of the details of the bandstructure and may be suited 
to cases where magnetic order is suppressed - like on a 
triangular lattice. Nevertheless, inter-site correlations, 
namely spin-spin interaction, can not be treated within 
the present approach. This is of particular interest in 
connection with superconductivity arising in doped Mott 
insulators. A study of superconductivity is presented 
elsewhere. 
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APPENDIX A: AVOIDED LEVEL CROSSING IN 
TWO MINIMAL MODELS 

If we consider a finite system in the z-direction com- 
posed of only L = 3 or L = 4 layers (but keeping the 
symmetry z <-> — z) the effective Schrodinger equation 
can easily be solved analytically. 

L = 3: The central layer may be almost localized with 
a strongly renormalized dispersion eo(e) = ZqE+Xo where 
Zq C 1. On the other hand, the two outer layers may be 
only weakly renormalized e-i(e) = z 2 e + \i = £i(e) with 



by V 



1. The hybridization (interlayer hopping) is given 
= —tzgzi. This leads to the following bands 



Ei = (e + ei)/2- ^ - e ) 2 /4 + 
E 2 = ei, 



E 3 = (eo + ei)/2 + ^(e 1 -e ) 2 /4 + 2V 2 . 

The renormalization of E\ and E3 is strongly e- 
dependent. Interestingly, one band is unchanged by the 
finite hybridization. The corresponding cigenfunction 
-02 = (1)0, — 1)/\/2 has no weight in the center. Fur- 
thermore, in the e-range where E\^ « E 2 , tpi,3 nas also 
only little weight in the center. The type of avoided level 
crossing in this minimal model corresponds to the type 
observed in heterostructures with even N (see Fig. [S]). 
L = 4: In this case the two central layers are character- 

£1/2 (e) and the outer 
= 63 /2(e). The inter- 
layer hopping matrix elements are given by Vq = —tz 2 ^ 2 



ized by e_x/ 2 (e) = zf /2 e + Ai /2 
layers by e_ 3 / 2 (e) 



/2 

Z 3/2 £ 



A 



3/2 



and V\ = — tei/223/2 which leads to the following bands: 
Vo 



Ei 



E 2 = e-^- 



Es 



Ei 




Ae + F 0l , y2 



Ae-Fn 



Here we use e = (ei/2 + £3/2)/2 and Ae = £1/2 — £3/2- The 
e-dependence of these bands is shown in Fig. [S] The type 
of avoided level crossing is generic for heterostructures 
with odd N. 



APPENDIX B: DRUDE WEIGHT AND 
TWISTED BOUNDARY CONDITIONS 

The Drude weight can be obtained from the ground- 
state energy of a system with twisted boundary condi- 
tions. The idea originally formulated by Kohn^L is to 
consider periodic boundary conditions in e.g. a;-direction 
and thread the system with a flux $. The flux can be 
represented by a vector potential A = (hc/e)Qx./N x (lat- 
tice constant a = 1) where N x is the number of lattice 
sites in x-dircction. We choose units such that h = c = 1 . 
Then, due to the Peierls phase, the hopping in a;-direction 
obtains an additional phase factor t — > t exp(±j$/A 2; ) 
which changes the free dispersion 

£fe ($) = -2t [cos(k x + + cos(k y )} . (Bl) 

The Drude weight can be calculated from the ground- 
state energy density E G (&) as^L^ 2 ^ 

d 2 E G (<S>) 



D = ixe Nt 



(B2) 

$=0 

has to be taken before 
00. In order to 



We note here that the limit $ 
going to the thermodynamic limit N : 
estimate the dependence of Eg on the twist $ we use the 
SBMF approximation. Because the explicit dependence 
of the saddle-point equations [Eq. ([51 El [Hi])] on ^ i s m 
second order (the ground state does not carry any cur- 
rent) we find that the linear change of the mean-fields in 
$ vanish 



ddi 




dni 




d\i 






$=0 


~d® 






$=0 



0. 



(B3) 



The calculation of D is thus considerably simplified be- 
cause only the explicit dependence on $ enters 

8 2 e g (<p) 



D 



2 

ire 



«9$2 



*=0 



E 

hif 



erf 



dZ, 



Ve k \ 2 f {E k „) 



(B4) 
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where n£° h is given in Eq. ([24]) and Z kv in Eq. (19]). This 
expression is equivalent to Eq. (|30[) . Furthermore, in sec- 
ond order perturbation theory we find that 



dZk v 



d 2 E kL , _ 
del ^ 



Ei. 



(B5) 



which shows that the nonuniform rcnormalization of the 
quasiparticle dispersion is directly connected to the spa- 
tial variation of zi . 

The integration necessary for the second term in 
Eq. (|B4[) is most conveniently done by introducing the 



weighted density of states^ 



(B6) 



which for a square lattice can be expressed in terms of 
complete elliptic integrals of the first and secon kind 



N v (e) 



St 



El 



,4tJ J \4t 
for lei < At and zero otherwise. 
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